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1. Introduction 

One of the greatest challenges in QCD is the calculation of higher orders in perturba- 
tion theory. Phenomenologically, these are important because the relatively large value of 
a s at accessible energies implies that many orders of perturbation theory are required in 
order to make precise quantitative tests. Theoretically, one expects the coefficients of the 
perturbative series for many QCD quantities to diverge factorially, and the rates of these 
divergences may cast light on issues in nonperturbative QCD, such as the existence and 
magnitudes of condensates and higher-twist effects [1] 

On the other hand, whilst progress in the exact calculations of higher-order terms in 
perturbative QCD series has been startling, with many new multi-loop results having re- 
cently become available [2] , existing perturbative techniques may not enable much further 
progress in exact calculations to be made in the near future. Thus various approximate 
techniques and numerical estimates may have a useful role to play. Among these, one may 
mention exact calculations of certain perturbative coefficients in the large- A^ limit, and 
the emerging lore of renormalons [1]. Also of potential use in QCD are Pade Approximants 
(PA's), as described in section 2 of this paper, which have previously demonstrated their 
utility in applications to problems in condensed-matter physics and statistical mechan- 
ics [3]. In recent years, these have been applied to obtain successful numerical predictions 
in various quantum field theories, including QCD, and justifications for some of these 
successes have been found in some mathematical theorems [4] on the convergence and 
renormalization-scale invariance of PA's. These theorems apply, in particular, to pertur- 
bative QCD series dominated by renormalon singularities, and in the large-,9o limit. 

Based on these theorems, a new method was introduced [5] for estimating the next- 
order coefficients in perturbative quantum field theory series on the basis of the known 
lower-order results and plausible conjectures on the likely high-order behaviour of the series, 
as also reviewed in section 2. This method "corrects" the conventional Pade Approximant 
Prediction (PAP) of the next term in the series by using an asymptotic error formula, 
providing improved predictions that we call Asymptotic Pade Approximant Predictions 
(APAP's). 

APAP's have already provided successful predictions for the perturbative coefficients 
in the subsequent calculation of the four- loop (3 function in QCD, as discussed in section 
3, and have also provided interesting results in N = 1 supersymmetric QCD (SQCD) [6]. 
The purpose of this paper is to provide a more complete account of these predictions, to 
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show how their accuracy may be improved in certain cases by a judicious weighting over 
negative numbers of flavours Np, and to extend these predictions to five loops in QCD in 
sections 5 and to SQCD in section 6. We also discuss analogous predictions for the QCD 
anomalous quark mass dimension in section 7, where "regular" APAP gives very good 
results, but the new weighting method does not improve matters. In section 8 we consider 
Abelian gauge theories, with less successful results. 

Before deriving these predictions, there is a technical issue which should be clarified, 
that may also illuminate an interesting physics point. As a general rule, (3 functions 
are scheme-dependent beyond one loop, and a theory with a single perturbative coupling 
constant g, such as QCD, is scheme-dependent beyond two loops, if one considers analytic 
redefinitions of g. In particular, the QCD (3 function can be transformed to zero beyond 
two loops, by making a suitable choice of renormalization scheme 2 . In our analysis of 
QCD, we use the MS scheme, and in iV = 1 SQCD we favour the DRED scheme 3 . 
The successes of the APAP procedure indicate that asymptotia and the convergence of 
PAP's are remarkably precocious in these schemes. In the SQCD case, there exists an 
alternative scheme (NSVZ) [7], associated with the Wilsonian action, in which there is 
an all-orders relation between f3 g and the quark anomalous dimension The NSVZ 
scheme differs perturbatively from DRED [8], and therefore provides a distinct test for 
the APAP method. We compare predictions for f3 g in both DRED and NSVZ, finding 
that they are less compelling in the latter case: perhaps minimal subtraction schemes are 
more amenable to Pade techniques? If so, it would be interesting to fathom the reason. 
As already noted, these techniques are not so successful for the quark mass anomalous 
dimension, or for Abelian theories. Perhaps these instances also provide clues when and 
why the Pade magic works. 

2. Formalism 

We start by recalling relevant aspects of the formalism for PA's and APAP's, and 
establishing our notation. For a generic perturbative series 



In fact, it can even be transformed to zero beyond one loop by a non- analytic redefinition 
of g involving In g: such redefinitions are associated with the Wilsonian action in supersymmetric 
theories. 

We recall that DRED corresponds to minimal subtraction in conjunction with regularisation 
by dimensional reduction. 




(2.1) 



n=0 
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the Pade approximant [N/M](x) is given by [3] 



ao + a\x + ■ ■ ■ cinx n 



b + bxx H b M x 

with b = 1, and the other coefficients chosen so that 

[N/M] = S + 0(x N+M+1 ). (2.3) 

The coefficient of the x N+M+1 term in (2.3) is the PAP estimate S$f M+1 of S N+M +i- If 
the perturbative coefficients S n diverge as n! for large n, it is possible to show [4] that the 
relative error 

cPAP c 
«->jy+M +l ~ ^N+ M+l 

'7V+M+1 

has the asymptotic form 



r _ "N+M + l -iv+iw+i , , 

On+m+i = o l 2 - 4 J 



M\A M 
Li 



Sn+m+i ^ "Ti (2.5) 



[N/M] 

as — > oo, for fixed M, where 

^[jv/m] = N + M + aM + b (2.6) 

and .A, a, 6 are constants. This theorem not only guarantees the convergence of the PAP's, 
but also specifies the asymptotic form of the corrections. 

The idea of APAP's is to fit the magnitude of this asymptotic correction using the 
known low-order perturbative coefficients, and apply the resulting numerical correction to 
the naive PAP's. In the applications discussed in this paper, we work with [0/1], [1/1] and 
[2/1] PA's, so that M = 1 throughout. For example, four-loop predictions are obtained as 
follows. In the case N max = 2, the [1/1] Pade leads to the naive PAP ^ Ap = S$/Si. The 
improved APAP estimate is then given by 

cPAP 

s 3 APAP = ^ 

1 + 03, 

where, motivated by its appropriateness in </> 4 field theory, we choose a + b = in the 
QCD application discussed in the next section, and A is then determined by comparing 
S2 to S2 AP = Si /So- Alternatively, we could have chosen a value of A and determined 
a + b from 62- However, as we shall see, when we go to five loops, knowledge of 82 and £3 
enables us to fit both A and a + b simultaneously. 
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3. Application to the Four-Loop (3 Function in QCD 

The APAP method was applied in [5] to estimate the four-loop QCD (3 function 
coefficient f3 3 , on the basis of the lower-order terms 



Po = ^-Ca — ^T F N F , 



Pi = fC A - AC pTpN p - fC A T F N F , 

a 2857 n% i at \on 2r r 205 n r> t 1415 ^2^ 

P2 — -54-^ A + ^pJ-F q-^F^A-I-F 27~ L/ A 1 F 



(3-1) 
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known before the appearance of the explicit four-loop calculation [9]. The quadratic 
Casimir coefficients C a and Cf for the adjoint and fundamental representations are given 
for the case of SU(Nc) by 

f 



C A = N c , 



N 2 

Of — 



2N C 



(3.2) 



and we assume the standard normalisation so that T F = ^ . We denote by Na the number 
of group generators, so that for SU(Nc) we have Na = N c — 1. 

We recall that j3 3 is a polynomial in the number of flavours N F : 



(3 3 = A 3 + B 3 N F + C 3 N 2 F + D 3 N F , 



(3.3) 



where D 3 = 1 .499 (for N c = 3) was already known from large-A F calculations. To justify 
applying the estimate (2.5), we assume that the f3 n ~ n\ for large n, as discussed in [5]. 
The predictions for A 3 , B 3 , C 3 resulting from fitting the APAP results for < N F < 4 to 
a polynomial of the form (2.5) are compared to the exact results in the first columns of 
Table I. 

The exact four-loop coefficient of the QCD f3 function for Nc colours is taken from 
the calculation of [9] , which was published after the APAP estimate: 
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where (3 = ((3) = 1.2020569- • •. The quartic Casimir coefficients in (3.4) are given for 
SU(N C ) by 



_N*{N* -l)(AT2 +36) 

a A a A — ^ 7 a F a A — ■ 

rl ab c d rl ab cd _ (N%-1)(N£.-6N*+18) 



(3.5) 



For Nc = 3 one obtains 

p 3 ~ 29243.0 - 6946.30N F + 405.089Nf + 1.49931-Nf., (3.6) 

whereas (3% is given by the coefficients shown in Table I when one omits the quartic Casimir 
contributions. 

These quartic Casimir terms appear for the first time at four-loop order. They are 
analogous to the light-by-light scattering terms in (g— 2) M , and PA-based techniques cannot 
estimate them on the basis of lower-order terms with different group-theoretical factors. 
Such terms are known to be important in (g — 2) M , but were relatively unimportant in 
previous perturbative QCD applications. In the case of /3 3 , they turn out to be about 15 
to 20 % for small N F , but are non- negligible for N F ~ 5. Setting these terms aside, the 
agreement between the predictions of [5] and the exact results of [9] is remarkable. The 
predictions we present in the rest of this paper should all be understood as applying to 
perturbative coefficients without the higher-order analogues of such quartic Casimir terms. 





APAP 


EXACT 


% DIFF 


WAPAP 


% DIFF 




23,600(900) 


24, 633 


-4.20(3.70) 


24, 606 


-0.11 


B 3 


-6,400(200) 


-6, 375 


-0.39(3.14) 


-6, 374 


-0.02 


c 3 


350(70) 


398.5 


-12.2(17.6) 


402.5 


-1.00 


D 3 


input 


1.499 




input 





Table I: Exact four-loop results for the QCD (3 function, compared with 
the original APAP's in the first column, and improved APAP's obtained 
from a weighted average over negative N F (WAPAP), as discussed in 
the text. The numbers in parenthesis are the error estimates from [5]. 
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Following [5], the same APAP method was applied in [6] to estimate the four- loop 
(3 function in SQCD. The agreement with known results was again encouraging, and the 
APAP provided a prediction a 2.4 for the unknown constant [8] in the four-loop SQCD 
P function, as discussed also in section 5. 



4. Weighted APAP's in QCD 

Before going on to make new predictions for QCD and SQCD at the five-loop level, we 
first draw attention to a refinement that offers an improvement on APAP's in the four-loop 
QCD case. As can be seen in Table I, the signs of the coefficients A3, B 3 , C 3 alternate. A 
corollary of this is that the APAP predictions for Np ~ 5 are sensitive to cancellations and 
relatively inaccurate. Conversely, the numerical analysis is relatively stable for (fictitious) 
N F < 0. We have observed empirically that more accurate predictions for the coefficients 
A3, B 3 and C 3 are obtained if one makes polynomial fits for some range of negative values 
of N F . Is there some systematic procedure that exploits this observation? The following 
is one method we have explored. 

We choose a range —N™ ax < N F < o over which we fit values of A using the APAP 
formulae of the previous section, and we determine the arithmetic mean of the correspond- 
ing values of A. We use this mean value of A to estimate (3 3 for each of the chosen values 
of N F , and fit to the polynomial form (3.3). We hypothesize that the most accurate results 
for the coefficients A3, B 3 , C 3 may be obtained when they contribute with equal weights to 
the fit: certainly, one cannot expect that any coefficient that has a small weight in the fit 
will be estimated reliably. For given iV™ 1 , the overall weights in the fit are A3, B 3 N F iax /2 
and C 3 N™ ax (2N™ ax + l)/6. We then estimate B 3 as follows. We take the two values of 
B 3 corresponding to the values of Np- ax for which the A3 and B3 weights are most nearly 
equal. Let us call these values of B3, B^ and B^ 2 \ and the corresponding weights B^^ 
and B^ ■ Our prediction for B 3 is then 

_ A 2 B^ + A 1 B™ 
B " ~ A 1+ A 2 ' (ZL1) 

where Ai 2 = \B™ — A^ \- We estimate C 3 in similar fashion. Both the S3 
calculation and the C3 calculation yield a result for A3, obtained as in (4.1) : we take as 
our prediction for A3 the mean of these two values. 
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Table I shows in the column labelled WAPAP the results we obtain using this proce- 
dure. We see that the latter are significantly more accurate than the ones obtained using 
the APAP's for < N F < 4. The values of N^ ax selected by WAPAP are 7, 8 for B 3 and 
13, 14 for C 3 . 

Table II compares the WAPAP predictions obtained in this way with the known exact 
results (omitting quartic Casimir contributions) in QCD for various values of Nc- The 
agreement is certainly impressive, even compared with the APAP results shown in Table 
I. Since the numerical value of the coefficient C 3 is relatively small, corresponding (in the 
case Nc = 3) to the relatively large value Np- ax = 14 mentioned above, it is perhaps not 
surprising that the percentage error in the estimate of this coefficient is larger than for 
either A 3 or S3. 

Figure 1 displays graphically our resulting predictions for /3 3 , as a function of N F 
for the most interesting case Nc = 3. We plot the percentage relative errors obtained 
using various APAP-based estimation schemes: naive APAP's fitted with positive Np < 4 
(diamonds), naive APAP's fitted with negative Np > —4, WAPAP's compared to the 
exact value of /3 3 including quartic Casimir terms, and WAPAP's compared to /3 3 without 
quartic Casimir terms (crosses). We see that the latter are the most accurate for f3 3 in 
QCD. In Figure 2 we show the error in the WAPAP prediction for f3 3 as a function of N F , 
for Nc =3, 4, 5, 6, 7 and 10, once again omitting quartic Casimir terms from the exact 
result. The accuracy of these predictions is our best evidence for believing in the utility of 
the WAPAP method. 

To anticipate the obvious question: we have explored whether this WAPAP procedure 
gives significantly better results than the conventional APAP's for the other perturbative 
series considered in this paper, namely the SQCD (3 function and the anomalous dimension 
of the quark mass. As we discuss in sections 7 and 8, the remarkable success of the method 
at four loops is not repeated for other cases, but there is distinct evidence (provided by 
large- A^-expansion results) that WAPAP leads to more reliable predictions at five loops. 
However, we feel that the results in Tables I and II already provide ample motivation for 
the QCD WAPAP calculation of (3± described in the next section. 
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WAPAP 



exact 



% error 



N c = 2 



A 3 


4.88 x 10 3 


4,866 


0.42 


B 3 


-1.86 x 10 3 


-1,854 


0.48 


c 3 


174 


170.5 


2.0 


N c = 3 


A 3 


2.467 x 10 4 


24, 633 


0.13 


Bs 


-6.383 x 10 3 


-6, 375 


0.13 


Cs 


405 


398.5 


1.6 


Nc = 4 


As 


7.790 x 10 4 


77, 852 


0.06 


B 3 


-1.521 x 10 4 


-15,210 


0.03 


Cs 


729 


717.2 


1.6 


N c = 5 


As 


1.901 x 10 5 


190,068 


0.04 


Bs 


-2.976 x 10 4 


-29,800 


-0.12 


C 3 


1.14 x 10 3 


1,127 


1.6 


N c = 6 


As 


3.943 x 10 5 


394, 125 


0.03 


Bs 


-5.149 x 10 4 


-51,580 


-0.17 


Cs 


1.65 x 10 3 


1,627.5 


1.6 


N c = 10 


As 


3.043 x 10 6 


3,041,089 


0.05 


Bs 


-2.388 x 10 5 


-239, 384 


-0.25 


Cs 


4.62 x 10 3 


4,540 


1.7 



Table II: Comparison of WAPAP and exact results for the exact 4-loop 
(3 function in QCD (omitting quartic Casimir terms), for various values 
ofN c . 
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w WAPAP 

vs |8 3 (N F ) w. quartic Casimirs 



Fig. 1: Predictions for (3$, as function of Np, for Nc = 3. The percentage 
relative errors are obtained using various APAP-based estimation schemes: 
naive APAP's fitted with positive Np < 4 (diamonds), naive APAP's fitted 
with negative Np > —4, WAPAP 's compared to the exact value of [3 3 including 
quartic Casimir terms, and WAPAP 's compared to (3s without quartic Casimir 
terms (crosses). 
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Fig. 2: The percentage relative errors in the WAPAP prediction for (3$ (com- 
pared to the exact result with quartic Casimir terras omitted), plotted vs. Np 
for N c =3, 4, 5, 6, 7, 10. 



5. Five-Loop Predictions in QCD 

We now outline the application of the APAP method to estimate the five-loop (3 
function coefficients (3^ in QCD, using our knowledge of the corresponding /3 to (3^. The 
standard [2, 1] Pade leads to the estimate 

/3f AP = f • (5-1) 
10 



This is then corrected in a similar fashion to Eq. (2.7): 



oPAP 

/3 4 APAP = (5.2) 
1 + 04 



where, according to Eqs. (2.5), (2.6), 5 4 is given asymptotically by 

A A 

5 4 = = ^— -. (5.3) 

L [2/ i] 3 + a + 6 v ; 

To estimate £4 we therefore need to know both .4. and a + b. These can be deduced from the 
lower-order relative errors 62 and £3, as defined in (2.4) for which we use the asymptotic 
estimates (2.5): 

^ = -(l + a + 6), ^ = _(2 + a + 6) (5.4) 

from which we obtain A and a + b 4 . 

We now calculate the WAPAP for the five-loop QCD (3 function, which we parametrise 

as 

p 4 = A 4 + B 4 N F + C 4 N 2 F + D 4 Nl + E 4 Np. (5.5) 

Once again we can input the coefficient of the highest power in Np, which is given in this 
case by [10]: 

E 4 = -4T£ [(288C(3) + 214)C F + (480C(3) - 229)C A ] /243 . (5.6) 
using which we obtain the five-loop results shown in Table III. 



4 The fitted value of a + b is not necessarily close to the value zero assumed in the estimate of 
03 m QCD. 
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N c 


2 


3 


4 


5 


10 


A 4 (w. Q) 


1.48 x 10 5 


7.59 x 10 5 


2.77 x 10 6 


7.92 x 10 6 


2.31 x 10 8 


A 4 (w/o Q) 


6.41 x 10 4 


4.88 x 10 5 


2.06 x 10 6 


6.28 x 10 6 


2.01 x 10 8 


B 4 (w. Q) 


-5.51 x 10 4 


-2.19 x 10 5 


-6.39 x 10 5 


-1.50 x 10 6 


-2.28 x 10 7 


B 4 (w/o Q) 


-3.04 x 10 4 


-1.56 x 10 5 


-4.97 x 10 5 


-1.22 x 10 6 


-1.95 x 10 7 


C 4 (w. Q) 


6.96 x 10 3 


2.05 x 10 4 


4.68 x 10 4 


9.00 x 10 4 


7.07 x 10 5 


C 4 (w/o Q) 


4.69 x 10 3 


1.64 x 10 4 


3.93 x 10 4 


7.72 x 10 4 


6.23 x 10 5 


L> 4 (w. Q) 


-21.8 


-49.8 


-89.8 


-142 


-575 


At(w/o Q) 


-28.3 


-60.5 


-105 


-163 


-640 


S 4 (input) 


-1.15 


-1.84 


-2.51 


-3.17 


-6.43 



Table III: WAPAP's for the five-loop QCD (3 function, calculated both with (w. Q) and 
without (w/o Q) the four-loop quartic Casimir terms in (3^. The values of N"p ax used 
range between 5 and 117 in the w. Q case, and between 4 and 108 in the w/o Q case, being 
largest for large Nc and for D 4 . 

Notice that in Table III we include results corresponding to both the inclusion (w. 
Q) and the omission (w/o Q) of the quartic Casimir contributions to the four-loop coeffi- 
cients, obtained from (3.4). The former (latter) results should of course be compared with 
contributions including (excluding) such terms at five loops when (and if) such results be- 
come available. Of course, at five-loop order we may expect to encounter new higher-order 
Casimir terms, which should in any event be omitted in the comparison. We can only 
hope that such contributions are relatively unimportant, which is the case for the quartic 
terms in /3s for small Np- We anticipate that the percentage errors of the w/o Q estimates 
of the non-quartic terms in the coefficients are likely to be the smallest, whereas the best 
estimate of the full coefficients may be provided by the w. Q estimates. 

We show below the results obtained if we choose not to input the value of i? 4 , but 
rather predict that as well. As can be seen, the results for A 4 , B4 and C 4 , in particular, 
are very stable. Moreover, the prediction for E4 is encouragingly close to the true value, 
considering the extreme smallness of E4 compared to A4. 
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N c 


2 


3 


4 


5 


10 


A 4 (w. Q) 


1.45 x 10 5 


7.51 x 10 5 


2.75 x 10 6 


7.87 x 10 6 


2.30 x 10 8 


A 4 (w/o Q) 


6.38 x 10 4 


4.85 x 10 5 


2.05 x 10 6 


6.24 x 10 6 


2.00 x 10 8 


B 4 (w. Q) 


-5.53 x 10 4 


-2.20 x 10 5 


-6.41 x 10 5 


-1.51 x 10 6 


-2.29 x 10 7 


B 4 (w/o Q) 


-3.05 x 10 4 


-1.57 x 10 5 


-4.99 x 10 5 


-1.22 x 10 6 


-1.96 x 10 7 


C 4 (w. Q) 


6.72 x 10 3 


1.97 x 10 4 


4.50 x 10 4 


8.66 x 10 4 


6.81 x 10 5 


C 4 (w/o Q) 


4.52 x 10 3 


1.58 x 10 4 


3.79 x 10 4 


7.43 x 10 4 


5.99 x 10 5 


£> 4 ( w - Q) 


-28.3 


-93.8 


-226 


-389 


-1,730 


D 4 (w/o Q) 


-72.7 


-163 


-287 


-446 


-1,750 


£ 4 (w. Q) 


-0.974 


-2.03 


-3.07 


-4.06 


-8.73 


£4 (w/o Q) 


-1.61 


-2.56 


-3.45 


-4.33 


-8.64 



Table IV: WAPAP's for the five-loop QCD (3 function, calculated with 
and without the four-loop quartic Casimir terms, but without inputting 
the known exact values of E 4 . It is encouraging to compare the output 
values with the last row in Table III. The values of N^ ax used range 
between 5 and 81 in the w. Q case, and between 4 and 104 in the w/o Q 
case. 

It is not possible to state precise errors for the type of prediction discussed in this 
paper. We gave in [5] certain estimates of the uncertainties, which turned out to be in 
the right ballpark if quartic Casimir terms are omitted in the comparison, as reported 
in Table I. The appearance of such new quartic terms is characteristic of the type of 
theoretical 'systematic error' that cannot be foreseen. In the case of our f3 4 predictions 
in QCD, we draw the reader's attention to the differences between the w. Q and w/o Q 
entries in Table III, and to the differences between these and the corresponding entries in 
Table IV, obtained without using the known values of E 4 as inputs. The most accurate 
estimates of the full coefficients are likely to be the w. Q entries in Table III, but the 
uncertainties are unlikely to be smaller than these differences. 

6. Five-Loop Predictions in N = 1 Supersymmetric QCD 

We begin with the SQCD f3 function in the DRED regularisation scheme, where the 
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first four coefficients are given by [8] : 



(3 = 3N C - N F , 



Pi = QN C - 
= 21N C 



4iV c 



2im 



N c 



N F 



-9 



N F - 



N c 



-AN C 



Ps = A 3 + B 3 N F + C 3 N 2 F + D 3 N F 

where Nc is the number of colours, and 
A 3 = (6 + 36a)N c 

B 3 = -36(1 + a)N c + (34 + 12a)N c + + 



(6.1a) 
(6.16) 
(6.1c) 

(6.1d) 



N c N c 

,'62 o \ , , 100 6k -20 



(6.2) 



3iV c 



Here k = 6C3 and a is a constant which has not yet been calculated exactly. Notice that 
there are no quartic Casimir contributions in the SQCD case 5 . The APAP method was 
used in an earlier paper [6] to obtain the estimate a ~ 2.4. 
Proceeding now to five loops, we write 



p 4 = A 4 + B 4 N F + C 4 N F + D 4 N F + E 4 N F . 



(6.3) 



As in the QCD case we can input the true value of E 4 provided by a recent large N F 
calculation [11], and given by 



E 4 = - [2N c Cs - (1 + 2Cs)/(2iV c )] . 



(6.4) 



We choose to calculate the WAPAP predictions both with and without this input. This 
also enables us to explore the sensitivity of the resulting prediction for E 4 to variations in 



5 Their absence may be understood as a consequence of the fact that the (5 function vanishes 
beyond one loop for an arbitrary N = 2 supersymmetric theory. We are unable, however, to 
comment on the possible appearance of quartic and higher-order Casimir terms at the five-loop 
level. 
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a. Assuming a = 2.4, we obtain the results shown in Table V, whereas the results with the 
known values of E 4 not input are shown in Table VI. The qualitative agreement between 
the predicted values of E 4 in the last row of Table VI and the exact values in Table V is 
good. We note that the WAPAP process is crucial for this agreement, in that the output 
E 4 is quite sensitive to the value of Np- ax used, which is fixed by the WAPAP criterion. 
We see that the output values of A4, B 4 , C 4 and D 4 are quite stable, which is perhaps 
to be expected in view of the small numerical values of E 4 . The differences between the 
results obtained with and without the input exact value of E 4 provide some indication of 
the uncertainty in the predictions. We expect, naturally, the case with input E 4 to be the 
more accurate. 



N c 


2 


3 


4 


5 


10 


A 4 


1.48 x 10 4 


1.13 x 10 5 


4.78 x 10 5 


1.46 x 10 6 


4.69 x 10 7 


B 4 


-1.05 x 10 4 


-5.85 x 10 4 


-1.91 x 10 5 


-4.72 x 10 5 


-7.70 x 10 6 


c 4 


3.25 x 10 3 


1.29 x 10 4 


3.21 x 10 4 


6.42 x 10 4 


5.29 x 10 5 


D 4 


-109 


-307 


-583 


-936 


-3.87 x 10 3 


E 4 (input) 


-3.96 


-6.64 


-9.19 


-11.7 


-23.9 



Table V: WAPAP's for the five-loop SQCD (3 function, assuming a = 2.4. The values of 
jyrnax use( ^ ran g e between 3 and 37. 



N c 


2 


3 


4 


5 


10 


A 4 


1.46 x 10 4 


1.12 x 10 5 


4.73 x 10 5 


1.45 x 10 6 


4.64 x 10 7 


B 4 


-1.04 x 10 4 


-5.87 x 10 4 


-1.91 x 10 5 


-4.74 x 10 5 


-7.73 x 10 6 


c 4 


3.16 x 10 3 


1.25 x 10 4 


3.11 x 10 4 


6.21 x 10 4 


5.12 x 10 5 


D 4 


-134 


-400 


-767 


-1.24 x 10 3 


-5.12 x 10 3 


E 4 


-2.44 


-4.53 


-6.33 


-8.03 


-16.1 



Table VP WAPAP's for the five-loop SQCD (3-function, again assuming a = 2 A, 
but without the exact values of E 4 as input. The values of ' N™ ax used range between 
4 and 61. 
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0.0 



N c = 3 




-10.0 



12.5 



-2 2 

a 

Fig. 3: The WAPAP result for E 4 plotted against a, for -3 < a < 3. 

The value a = 2.4 used above was itself based on an APAP calculation [6]. It behoves 
us, therefore, to explore the sensitivity of our results to the precise value of a. In Fig. 3 we 
plot the WAPAP result for E 4 against a, for —3 < a < 3. We see that for this range there 
are two values of a corresponding to E 4 = i?| xact , namely a —0.9 and a « 1.4. Given 
the fact that in general we would expect E 4 to be the least- well determined coefficient, we 
consider this result to be reasonably consistent with our previous prediction that a ~ 2.4. 
It should be noted that our predictions for A 4 ■ ■ ■ D 4 are also sensitive to the precise value 
of a. 

We turn now to the alternative NSVZ prescription for the SQCD /3-function, given by 
the following exact formula [7] which relates f3 g to the quark anomalous dimension, 7 9 : 



0\ 



■NSVZ 
9 



g- 



16tt 



N F - 3N C - 2N F 7, 



NSVZ 



1 -2N c g 2 (l6ir 2 y 1 



(6.5) 



Note the overall minus sign, in accordance with our conventions here. Using (6.5) and the 
result for 7^ svz given in [8], we obtain: 



A) = 3A^ C - N F , 



(6.6a) 



16 



Pi = 6N C - 
(3 2 = 12N C - 



AN C - 



N c 



N f 



12N, 



c 



N 2 C 



N F 



N c 



2N C 



Ni 



[3 3 = A 3 + B 3 N F + C 3 N 2 F + D 3 N F 



(6.66) 
(6.6c) 

(6.6c?) 



where 



A 3 = 24N C 

B 3 = -AON c + 3(W C - -f- + 



C 3 = (2k + 14) N c - 24 
2 



N c N c 
2k- 10 



D 3 = 2N C - 



N c 



(6.7) 



In this case there is no undetermined parameter a: we know [8] 7^ svz through three loops 
and hence /3^ svz through four loops. 

It is possible to argue [12] on the basis of the nature of the coupling-constant redefi- 
nition connecting the two schemes that -^red anc j -yNSVZ are ^ e same at leading order 
in N F . Hence, if as before we write 



p 4 = A 4 + B 4 N F + C 4 N F + D 4 N F + E 4 N F , 



(6.8) 



we can find E 4 , as we did in the DRED case, from the large- N F results in [11]. The result 
is: 

£ 4 = 2[l-2C(3)](iV c -l/iV c ). (6.9) 

We also have, as is evident from (6.5), that A 4 = 48N C , providing an additional check on 
our calculation 6 . Our WAPAP results are shown in the Tables VII and VIII, for the cases 
with and without E 4 input. Also shown in the second row of Table VIII are the exact 
results for A 4 . 

6 We could, of course, input both A4 and E 4 , but we choose instead to compare the WAPAP 
results for all the five-loop coefficients with the corresponding ones with E4 input. 
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N c 


2 


3 


4 


5 


10 


A 4 


1.68 x 10 3 


1.04 x 10 4 


4.44 x 10 4 


4.99 x 10 5 


4.42 x 10 6 


B 4 


-1.25 x 10 3 


-7.87 x 10 3 


-2.63 x 10 4 


-6.56 x 10 4 


-1.08 x 10 6 


c 4 


750 


3.11 x 10 3 


7.87 x 10 3 


1.58 x 10 4 


1.32 x 10 5 


D 4 


-6.0 


-90.1 


-163 


-516 


-938 


E 4 (input) 


-4.21 


-7.49 


-10.5 


-13.5 


-27.8 



Table VII: WAPAP's for the five-loop NSVZ f3 function, with the exact values of E 4 used 
as input. The values of N™ ax used range between 3 and 26. 



N c 


2 


3 


4 


5 


10 


A 4 


1.49 x 10 3 


1.05 x 10 4 


4.33 x 10 4 


1.42 x 10 5 


4.45 x 10 6 


A 4 (exact) 


1.536 x 10 3 


1.166 x 10 4 


4.915 x 10 4 


1.500 x 10 5 


4.800 x 10 6 


B 4 


-1.13 x 10 3 


-7.80 x 10 3 


-2.65 x 10 4 


-6.64 x 10 4 


-1.09 x 10 6 


c 4 


612 


2.87 x 10 3 


7.35 x 10 3 


1.48 x 10 4 


1.23 x 10 5 


D 4 


-75.2 


-241 


-462 


-742 


-3060 


E 4 


-13.0 


-13.3 


-15.8 


-27.9 


-50.6 



Table VIII: WAPAP's for the five-loop NSVZ (3- function, with E 4 not input. The values 
of Np- ax used range between 4 and 25. 

We see that the WAPAP's are in general in good agreement with the exact result for 
A 4 in the NSVZ scheme, at the 10% level. Although encouraging, these results are not 
quite as compelling as the ones for the DRED scheme. This is at first sight surprising, 
given the form of (6.5), which appears at first sight to be close to the rational function 
form of the PA's. However, as mentioned in the Introduction, perhaps minimal subtraction 
schemes are more amenable to Pade techniques. The anomalously poor result for A 4 in 
Table VII is caused by the fact that the error £4 is close to —1 in this case, for the N™ ax 
values corresponding to the determination of D 4 . The reason the result for D 4 is not also 
anomalously large is that the two values from which the weighted average is taken are both 
numerically large but with opposite signs. Thus we cannot rely on either the A 4 or the D 4 
prediction for Nc = 5. With this exception, A 4 comes out reasonably close to the exact 
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result. This means, of course that the predictions for £? 4 • • • D 4 will not change much if we 
input A 4 as well as E4. 

Analogously to the five-loop QCD case discussed in the previous section, we take 
the differences between the entries in Tables VII and VIII as lower limits on the possible 
uncertainties in our five-loop NSVZ predictions. 



7. The Quark Mass Anomalous Dimension in QCD 

We now consider the quark mass anomalous dimension 7 in QCD, defined as 



7 = 



dhim 



- = -70a - 71a 2 - 7 2 a 3 



73a 4 - 74 a 5 +0(a 6 ), 



(7.1) 



where a = a s /n. The 4-loop coefficient 73 was recently computed in [13] [14] and the full 
exact results for the coefficients 7™ for n = 0, 1, 2, 3 are given by 

70 = \ [3Cf] 

71 = J6 [f + ^-CfCa — ^-CfTfNf] 



72 64 



129^3 129 r* 2 n 1 11413/1 r* 2 
T°F —^F^A "T 108 



+ C 2 TpNp(-46 + 48C 3 ) + C F C A T F N F (-f£ - 48C 3 ) - ±§CpT 2 N 2 



73 



1 

256 



C£ (-T - 336C 3 ) + C%C A (i^ + 3I6C3) 



+ Cf^a (" 



34045 
36 



- 152C 3 + 440C 5 )+C F Ci( 



70055 
72 



+ 



1416 



Cs " 440C5) 



+ CpTpNp (-^P + 552C 3 - 480Cs) + C 2 F C A T F N F (-^ + 368Cs - 264C 4 + 8OC5) 



+ C f C 2 a T f Nf 



65459 
162 



2684 



C 3 + 264C 4 + 400C 5 ) 



+ C|T|iV| - I6OC3 + 96C 4 ) 
+ C F C A T 2 N 2 (±fp + I6OC3 - 96C 4 ) + C F T F N F (- 

JabcdJabcd r ]abcd r ]abcd 

+ , (-32 + 240C 3 ) + iV F F / F (64 - 480C 3 ) 



664 , 128 
81 9 



Q 



(7.2) 

where for SU(Nq) the quadratic and quartic Casimirs are as defined in (3.2), (3.5), and 
T F = \ as before. In addition, dQ is the dimension of the quark representation, so 
that d Q = N c for SU(N C ), and we have C4 = C(4) = 1.0823232 •• • and Cs = C(5) = 
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1.0369277- 



For Nc = 3, we have 



7o = 1 

72 = ± [1249 + _ 160 Ca) ^ _ ig,^] (7 3) 

1 f 4603055 j_ 135680 /- oonn/- _i_ f 91723 34192 /- , oon/ _L 18400 /- \ rvr 

73 = 256 [ 162 + ^7-^3 ~ 88OOC5 + ( 27 S - S3 + OOUQ + — g— C5J 

1 /5242 , 800 a 160a \ at2 , / 332 , 64a \ Ar3 1 
+ I "243" + — « - — C4j A/ F + (-243 + 27C3) ^fJ , 



which have the numerical values 



70 = 1 

71 « 4.20833 - 0.138889N F 

72 « 19.5156 - 2.28412iV> - 0.0270062Nf 

73 « 98.9434 - 19.1075iV> + 0.276163A^| + 0. 00579322 iV|. 



(7.4) 



Omitting the quartic Casimir contributions, one obtains 

73 = 96.4386 - 18.8292iV> + 0.276163iVf + 0. 00579322 iVf,. (7.5) 

and we shall now compare (7.4) and (7.5) with APAP's. 

It transpires that the WAPAP procedure does not work so well here. The most 
accurate results for both and C~ are obtained for small N^ ax . This is reasonably 
consistent with the WAPAP behaviour in the C~ case: here, the weight difference C] W — 
AJ W never changes sign, but is smallest at Np ax = 2 on the edge of the range. However, 
the WAPAP criterion for B% leads to values of Np ,ax which start at 9 for Nc = 2 and 
increase with Nc- Nevertheless, as in the previous sections, it seems sensible to match 
at negative N F , and spectacular results are obtained if we simply take N^ ax = 4 (with 
_pj-rnax < < throughout, as can be seen from Table IX, where numerical predictions 
for the coefficients in the parametrization 

73 = Al + BJN F + cqN% + DJN 3 F , (7.6) 

are given both without (w/o Q) and with (w. Q) quartic Casimir contributions. It should 
be noted that we have used as input the exact result for Dj, which is contained in (7.3). 
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N c 


2 


3 


4 


5 


20 



APAP 


16.1 


97.9 


328 


822 


2.18 x 10 5 


w/o Q 


15.4 


96.4 


327 


825 


2.23 x 10 5 


w. Q 


16.0 


98.9 


334 


840 


2.26 x 10 5 


B 1 


APAP 


-5.14 


-20.0 


-49.3 


-98.0 


-6.39 x 10 3 


w/o Q 


-4.70 


-18.8 


-47.1 


-94.2 


-6.27 x 10 3 


w. Q 


-4.77 


-19.1 


-48.0 


-96.2 


-6.43 x 10 3 


cq 


APAP 


0.065 


0.224 


0.478 


0.828 


17.5 


exact 


0.111 


0.276 


0.504 


0.796 


13.0 


DJ 


input 


3.26 x 10" 3 


5.79 x 10-3 


8.15 x 10-3 


0.0104 


0.0433 



Table IX: Four-loop quark mass anomalous dimension in QCD: APAP's 
for fixed jVj^ ax = 4 are compared with the exact values both without (w/o 
Q) and with (w. Q) the quartic Casimir terms. 

It can be seen that in all cases the APAP estimate is quite accurate over a wide range 
of Nc- In most cases, the APAP estimate is closer to the exact result without the quartic 
Casimir contribution (w/o Q), but in any case the quartic Casimir contribution to 73 is 
smaller than in the case of the QCD ft function. 

We now go on to discuss the five- loop APAP estimate of 7. We parametrize the 
five-loop quark mass anomalous dimension 74 in the form 

74 = A\ + BJN F + CJN 2 F + DJN 3 F + EjN F , (7.7) 

where the value of Ej can be derived from [15]: 

EJ = C F T| (-65/5184 - 5C(3)/324 + tt 4 /3240) . (7.8) 
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We use the full 73 as input, including the quartic Casimir contribution. As we argued in 
the case of the QCD f3 function, we expect our five-loop estimate to include the effects 
of contributions involving such quartic Casimir terms, but not the effect of new Casimir 
terms making a first appearance. Once again we choose N F %ax = 4 to derive the results 
shown in Table X. 



N c 


2 


3 


4 


5 


20 


A 4 (w. Q) 


56.0 


530 


2.41 x 10 3 


7.63 x 10 3 


8.37 x 10 6 


A 4 (w/o Q) 


50.5 


493 


2.27 x 10 3 


7.22 x 10 3 


7.97 x 10 6 


S 4 (w. Q) 


-23.3 


-143 


-483 


-1.22 x 10 3 


-3.33 x 10 5 


B 4 (w/o Q) 


-21.7 


-135 


-457 


-1.15 x 10 3 


-3.12 x 10 5 


C 4 (w. Q) 


1.70 


6.67 


16.8 


33.7 


2.29 x 10 3 


C 4 (w/o Q) 


1.64 


6.44 


16.0 


32.0 


2.14 x 10 3 


£> 4 (w. Q) 


8.12 x 10" 3 


0.037 


0.0891 


0.165 


4.31 


D 4 (w/o Q) 


8.88 x IO- 3 


0.037 


0.0831 


0.148 


3.48 


E4 (input) 


-4.80 x 10~ 5 


-8.54 x 10~ 5 


-1.2 x 10~ 4 


-1.54 x 10~ 4 


-6.39 x 10~ 4 



Table X: APAP's for the five-loop quark mass anomalous dimension in QCD, calculated 
with and without the four-loop quartic Casimir terms. 



8. Abelian Gauge Theories 

All of the previous sections have dealt with APAP predictions for non- Abelian theories. 
It is natural to ask whether similarly accurate results can be obtained for the Abelian 
case. We address this question in this section, choosing as our example the fermion mass 
anomalous dimension with Np charged fermions, where good results were found in the 
non- Abelian case, as we saw in the previous section. A supplementary reason for choosing 
this example is that the next-to -leading- N p result is available, as well as the leading one. 

The results for 71 • • • 73 in the abelian case follow from Eq. (7.2) by setting 

labcd labcd r labcd r labcd 

C F = T F = l,C A = r F f A = 0, f f F = 1, (8.1) 

cLq d Q 
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so that: 



70 = 0.75 (8.2a) 

71 « 0.09375 - 0.2083iV F (8.26) 

72 « 1.0078 + 0.18279iV F - 0.08102iV F (8.2c) 

73 « -2.1934 - 1.7207iV F - 0.30143iVf + 0.03476iV F . (8.2d) 

Omitting the quartic Casimir term, we would instead have 

73 « -2.1934 + 0.2831iV F - 0.30143iVf + 0.03476A/|. (8.3) 

We can see at once that the miraculous success of the previous APAP prediction for 73 
will not be reproduced here. For N F = 0, the sign of 73 differs from the sign of 72/71- 
Moreover, 71 has a zero, and hence 7 3 4PAP has a pole, for Np ~ 0.45. Hence, we cannot 
hope to reproduce 73 for small values of \Np\. For large \Np\ the sign of 73 is still wrong, 
so the method fails in this region also. 

One easily verifies that this pessimism is confirmed by the results, and things do not 
improve at five loops. Then, as well as Ej as given in (7.8) , it is possible to derive from 
[16] the result for D 4 : 

DJ = -C 3 + -C5 - — + - 0.0804. (8.4) 

4 96 s 6 s 288 41472 v 1 

We notice now, however, that 72 has zeros, and hence 7 APAP has poles, for Np = —2.6 
and Np = 4.8. Consequently, we may expect that the results will be rather sensitive to 
the range of Np, if we match in a region including the origin. On the other hand, for large 
Np we have 73/72 ~ — 0.014AT|,, whereas Ej —0.001, so we also cannot expect good 
results at increasing \Np\. 

We leave it to the reader to convince her (him) self that we cannot expect to extract 
reliable predictions for A±---C±. We also record that the QED and SQED gauge (3- 
functions yield similarly unattractive results. Evidently, Abelian theories are less amenable 
to the APAP approach, for some unknown reason. 
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9. Conclusions 

We have presented results obtained from our APAP method for the four-loop and 
five-loop QCD /3-function coefficients, for the five-loop SQCD /3-function coefficients, and 
for the four- and five-loop quark mass anomalous dimensions in QCD. Particularly in the 
case of the QCD /3-function, and to some extent also for SQCD, particularly in the DRED 
scheme, a modified procedure for extracting the predictions for the various coefficients of 
powers of N F (WAPAP) gave improved results. In general, the four-loop results agree very 
well with the known results, giving us confidence in our predictions of the five-loop terms. 

Our four-loop QCD /3-function predictions [5] were confirmed very rapidly by an exact 
calculation [9]. Unfortunately, in view of the current limitations on the technology of 
exact perturbative calculations in QCD and SQCD, it may be some time before our five- 
loop predictions can also be tested directly. It would therefore be interesting to find 
alternative techniques that could be confronted or combined with APAP's. One possible 
complementary technique may be that of the large-A^ expansion. Unfortunately, it is 
the leading term in N F which is least well determined by the APAP approach, which is 
related to the poor results obtained in the Abelian case. It would be very interesting if 
the large-A^ methods could be extended to next-to-leading terms in this expansion for 
the non-Abelian case, in which case more comparisons and cross-checks could be made. 
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